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ABSTRACT 

We present the first observational limits to the ultra-high energy (UHE) neutrino flux 
from the Galactic Center, and from Centaurus A which is the nearest active galactic 
nucleus (AGN). These results are based on our "Lunar UHE Neutrino Astrophysics 
using the Square Kilometer Array" (LUNASKA) project experiments at the Australia 
Telescope Compact Array (ATCA). We also derive limits for the previous experiments 
and compare these limits with expectations for acceleration and super-heavy dark 
matter models of the origin of UHE cosmic rays. 
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1 INTRODUCTION 

Arrival directions of the UHE cosmic rays (CR) detected by 
the Pierre Auger experiment above 5.6 x 10^^ eV have been 
found to be statistically corre lated with positions of nearby 
AGN (| Abraham et al.ll2007al '). and a few of the arrival di- 
rections appear to be clustered around Centaurus A, our 
nearest active galactic nucleus at a distance of ~3.7 Mpc. 
This has led to speculation that Centaurus A may be re- 
sponsible for some of the UHE CR. However, the flux is 
extremely low, and so the nature of the sources of UHE CR 
remains at present unknown. 

As well as observing UHE CR directly, an alternative 
means of exploring the origin of UHE CR is to search for 
UHE neutrinos. Cosmic rays of sufficient energy will interact 
(e.g. via pion photo-production) with photons of the 2.725 K 
cosmic microwave background (CMB) radiation, with the 
resulting energy-loss producing a cut-off in the spectrum 
at around ^ fO'^° eV fro m a distant source ((Grcison 196& 
IZatsepin fc Kuzmirj|l966l ). These same interactions produce 
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Ensel et al.ll2001). As well 



as these cosmogenic neutrinos, UHE neutrinos are also ex- 
pected to be produced by acceleration and super-heavy dark 
matter (SHDM) sources of UHE CR, and some informa- 
tion on the CR spectrum at the sources is imprinted on 
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the spectrum of cosmogenic neutrinos (|Protheroell2004l ). Of 
course, neutrinos are not deflected by magnetic fields, and 
so should point back t o w here they were pro duced. See 
IProtheroe fc Clavl (|2004l ) and lFalcke et al.l (|2004l ) for reviews 
of UHE CR production scenarios and radio techniques for 
high-energy cosmic ray and neutrino astrophysics. 



1.1 The Lunar Cherenkov Technique 

In our present work we use the luna r Cherenkov technique 
(|Dagkesamanskii fc Zheleznvkhll 19891 ) . in which the Moon is 
used as a UHE neutrino target and Earth-based radio tele- 
scopes are used to detect coherent radio Cherenkov emis- 
sion produced by neutrino-induced cascades in the lunar re- 
golith. A high-energy particle interacting in a dense medium 
will produce a cascade of secondary particles which devel- 
ops an excess negative charge by entrainment of electrons 
from the surrounding material and positron annihilation in 
flight. The charge excess is proportional to the number of 
particles in the electromagnetic cascade, which in turn is 
propo rtional to t he energy of the primary particle. I Askar'vanI 
l| 19621 ) (see also lAskar'vanlll965l ) first noted this effect and 
predicted the Cherenkov emission process in dense dielectric 
media to be coherent at radio frequencies where the wave- 
length is comparable to or larger than the dimensions of 
the s hower, and this effect has been confirmed experimen- 
tally (jSaltzberg et al.ll200ll ). At wavelengths comparable to 
the width of the shower, the coherent emission is in a nar- 
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Figure 1. Contours of effective area (km^) as a function of UHE 
neutrino arrival direction for 10^'^ eV neutrinos for limb-pointing 
configuration (May 2008). The '+' marks the position of peak 
effective area; the Moon is at the center 5)=(0,0); telescope 
pointing direction (rj, C)=(0.183° ,0.183°). 



row cone about the Cherenkov angle, while for wavelengths 
comparable to t he shower length the coherent emission is 
nearly isotropic ({Alvarez- Munizll2006l ). 

The lunar Cherenkov technique aims to utilize the outer 
layers of the Moon, nominally the regolith which is a sandy 
layer of ejecta covering the Moon to a depth of ~10 m, as a 
suitable medium to observe the Askaryan effect. Since the re- 
golith is comparatively transparent at radio frequencies, co- 
herent Cherenkov emission from cascades due to sufficiently 
high-energy neutrino interactions in the regolith should be 
detectable as nanosecond-scale pulses by Earth-based radio- 
telescopes. 



2 DESCRIPTION OF THE EXPERIMENT 

The aim of the LUNASKA project is to develop further 
the lunar Cherenkov technique for UHE neutrino astron- 
omy, and to influence the design of the Square Kilometer 
Array (http://www.skatelescope.org/) so that UHE neu- 
trino observations may be possible. Our experiment was 
carried out at the ATCA which is an aperture synthesis 
telescope located at latitude —30° near Narrabri, Australia. 
It consists of six identical 22 m dishes of which we have 
used three with baselines ranging from ~100 m to ~750 m. 
The ATCA has a half-power beam-width that matches the 
lunar disk at 1.4 GHz, and it provided us with 600 MHz 
(1.2-1.8 GHz) of bandwidth. Being an array the ATCA also 
provided both strong timing discrimination against terres- 
trial radio- frequency interference (RFI), and gave a large 
effecti ve area and high sensitivity while seeing the entire 
moon lEkers et al.ll2009l ). In order to perform a search for 
nanosecond-duration lunar Cherenkov pulses, we had to 



build specialized hardware, including analogue de-dispersion 
filters as such pulses suffer dispersion in the Earth's iono- 
sphere - our experiment is the first to coherently correct 
for this before detection in real-time. To detect and store 
candidate events in real time we used field-programmable 
gate array (FPGA) based analog-to-digital converters devel- 
oped by the Australia Telescope National Facility, each of 
which could digitize and perform simple logic on two data 
streams at a rate of 2.048 GHz. The signal was passed to 
both a running buffer of length 256 samples and a real-time 
trigger algorithm. We triggered independently, with a maxi- 
mum rate of 1040 Hz, at each antenna. On fulfilling the trig- 
ger conditions, the buffer was returned to the control room 
and recorded. We calibrated the system sensitivity using the 
thermal emission from the lunar disk, and the system clocks 
at each antenna using correlated emission from the quasar 
3C273 . Ful l details of the exp eriment are given by Jamgi 
(120091 ) and I James et all (|20ld ). 

The observations described here cover two observing 
periods, February and May 2008. The February 2008 ob- 
servations were tailored to target a broad (> 20°) re- 
gion of the sky near the Galactic Center, both a poten- 
tial accelerator of UHE CR and also a potential source of 
UHE CR, gamma-rays and neutrinos t hrough its dark mat- 
ter halo. Based on simulation re sults (|james &: Protherod 
'2009a|;|j ames fc Pro thcroc '2009bh , we pointed the antennas 
at the lunar center in February to achieve the greatest to- 
tal effective aperture and sensitivity to an isotropic or very 
broadly-distributed fiux. Our May 2008 observing period 
targeted Centaurus A only, and by pointing the telescopes 
at the portion of the lunar limb closest to Centaurus A we 
achieved maximum sensitivity for this source. We recorded a 
total of 98307 3-fold coincidences within 4-/is windows, the 
majority of which, based on statistical arguments, could not 
come from random noise. Fitting these to far-field sources, 
~ 60 appeared to come from the direction of the Moon. 
However, our 0.5 nanosecond timing resolution allowed us 
to show that these were all either near-field in origin, or 
were identifiable under visual inspection as long-duration, 
narrow-band RFI signals with little or no time-structure, 
occurring during short periods of intense RFI. Thus we ex- 
cluded all candidates as being of lunar cosmic-ray or neu- 
trino origin; cosmic rays interacting with lunar regolith with 
favourable geometry could also produce observable lunar 
Cherenkov emission. 

We calculated the eff ective areas as described by 
iJames fc Protherod (|2009bh . This is shown in Fig. [1] for 
10'^'' eV neutrinos and the limb-pointing configuration, and 
we see that the sensitivity has a characteristic 'kidney' shape 
peaking at ~15° away from the Moon along the line extend- 
ing from the Moon's center to the telescope pointing posi- 
tion on the lunar limb. For the center-pointing configuration 
(not shown), the sensitivity pattern forms an annulus, with 
peak exposure around 15°-20° degrees from the Moon. For 
both configurations, the sensitivity pattern broadens with 
increasing primary particle energy as the increased strength 
of the pulses produced allows the telescopes to be sensitive 
to a wider range of interaction geometries. Combining the 
instantaneous aperture, e.g. as shown in Fig. [1] with the 
known telescope-pointing positions on the Moon and the 
Moon's position itself at the time of observing, allows us 
to a calculate the exposure [A x t\{Ev; oi, 5) (effective area- 
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time product) as a function of celestial coordinates {a,S). 
The exposure from our LUNASKA ATCA observations to 
UHE neutrinos is shown in Fig. [2] for 10^'^ eV. The nominal 
declination range of the ANITA observations is also given, 
and in this range the exposure of ANITA (not shown) dom- 
inates. The concentration of exposure about Centaurus A, 
and the broad Galactic Centre region (nominally Sagittarius 
A), both of which are outside ANITA's sensitive declination 
range, is due to the targeting of these regions in our exper- 
iment by a careful choice of observing times. 



3 RESULTS 

The experiment with the greatest exposure to UHE neu- 
trinos in the energy range above 10^^ eV applicable here 
is the An tarctic Impulsive Transient Antenna Ex periment 
(ANITA) (|Barwick et al.ll2006l : [Gorham et al.ll2010l ) but this 
experiment is only sensitive to the declination range —10° < 
S <-l-15° (indicated in Fig. [2]). The IceCube experiment is 
primarily sensitive to the northern hemisphere, but has nev- 
ertheless been able to set limits up to 10^^ eV for C entau- 
rus A and the Galactic Center l|Abbasi et al.l 120091 '). The 
Pierre Auger observatory, although designed as a cosmic 
ray detector, is also sensitive to neutrinos, and a prelimi- 
nary limit has been published for neutrinos from Centau - 
rus A with energies up to 10^° eV (|Tiffenberg et al.ll2009l ). 
Other experiments for which we can obtain a neutrino flux 
limit from Centaurus A above 10^^ eV ar e the pioneering 
experiment at Parkes (H ankins et al.|[l996l ). the Goldstone 
Lunar UHE Neutrino Experiment (GLUE) (|Gorham et al.l 
I2OO4I '). the Nu Moon observat i ons u sing the Westerbork 
interferometer (|Scholten et al.l |2009|). and the Radio Ic e 
Cherenkov Experiment fRICE) (|Kravchenko et al.1 \20odt ). 
The directional dependence of the exposures of the original 
Parkes experiment and GLU E were calculated previously by 
IJames fc Protherod (|2009bh . 

The NuMoon lun ar Cherenkov observa tions at Wester- 
bork are reported by IScholten et al. While neither 
Centaurus A nor the Galactic Centre was explicitly targeted, 
lunar Cherenkov experiments operating at an order of mag- 
nitude lower frequency than ATCA, such as the Westerbork 
array, ar e expected to be mo re sensitive at energies above 
10^^ eV l|Scholten et al. ''2009Y We simulated the sensitivity 
to isotropic 10'^'' eV neutrinos for the NuMoon experimen- 



tal configuration of two fan beams on different sides of the 
Moon, each with 4 frequency sub-bands of width 65 MHz 
centred on 123, 137, 151 and 165 MHz covering 1/3 of the 
Moon. Adopting a 10 m thick regolith with a denser sub- 
re golith layer, and using the triggering criteria as described 
by IScholten et al.l (|2009l ) we found a sensitivity a f a ctor o f 
^ 10 lower than that calculated by IScholten et ai] (jiooi). 
We are able to closely reproduce the published NuMoon re- 
sult by instead simulating one continuous bandwidth with a 
240 kjy threshold (artificially improving the sensitivity by 
a factor of 5) and letting the depth of the top regolith layer 

become infinite (gaining a furt her factor of 2). 

Using the method of iJames fc Protherod ()2009bl ) , but 
taking account of the fan beams of the Westerbork antennas, 
the frequency sub-bands and triggering criteria, our simu- 
lations show that for NuMoon the instantaneous neutrino 
directional aperture rises rapidly from zero at 0° from the 



Moon to a peak at ~ 40° away from the Moon, and then 
drops to zero at ~ 100° away from the Moon. To calculate 
the li mit for Centau rus A we obtained the observation dates 
from iBuitinkI l|2009i) . and on two of these dates the Moon 
happened to be ~ 45° from Centaurus A. Because only the 
observation dates (not the times) were available we assumed 
uniform lunar coverage by the two beams, but reduced by 
1/3 since each beam covered ~ 1/3 of the Moon. We calcu- 
lated the neutrino sky coverage of NuMoon by integrating 
the instantaneous neutrino directional aperture for the posi- 
tions of the Moon on the observation dates and obtained the 
NuMoon sensitivity for 10^^ eV neutrinos from Centaurus 
A. We obtained the sensitivity for the Galactic Centre using 
the same method. 

RICE was a Cherenkov radio experiment embedded 
in Antarctic ice at the South Pole and had a much lower 
neutrino energy threshold than the three lunar Cherenkov 
experiments, albeit with only a slowly-increasing exposure 
with neutrino energy. However, it had a very long observa- 
tion time of several years as compared to several days for 
the lunar Cherenkov experiments and this compensates for 
its lower instantaneous effective aperture. RICE was mostly 
sensitive to down-going neutrino events, and so to declina- 
tions 5 <0°. T he directional sensitivity pattern at lO'^^ eV 
( Bos sonll2008l ) shows that its sensitivity to Centaurus A and 
the Galactic Centre are respectively 1.0 and 1.4 times its 
average sensitivity to a source in the southern hemisphere. 
Taking this directional sensitivity as representative for all 
energies, and the energy dependence of th e expo sure to a 
diffuse flux as given by iKravchenko et ahl l|2006l ) we can 
thus calculate the exposure of RICE to these point sources. 
Since the isotropic exp osure is only given to lO^'^ eV in 
IKravchenko et all (|2006l ). in calculating the directional ex- 
posure we use a log-linear scaling and assume an isotropic 
exposure of 6.3 x 10^^ cm'^ s sr to 10^^ eV neutrinos. 

The individual and total exposures of the GLUE, RICE, 
NuMoon and LUNASKA experiments to UHE f from Cen- 
taurus A and Sgr A are given in Table [1] for 10^\ 10^^, 
and 10^"^ eV. Model-independent 90% confidence limits to 
the neutrino fiux at energy Ei, from putative point sources 
at positions (a, 5) can be obtained from the exposure us- 
ing E^F^{eJ) s; 2.3/[A X t\{E^; a, 5). The errors in the LU- 
NASKA exposures given in Table [T]refiect the in-quadrature 
addition of two dominant factors - uncertainties in the ab- 
solute detector thresholds, and the unknown properties of 
the lunar regolith. An esti mate of the forme r is described 
in our earlier publication (|James et al.ll2010l ). To estimate 
uncertainties due to regolith properties, we use the mea- 
sured variation in moon roc k radio absorption and density 
(|01hoeft fc StrangwavllT975h . The mare have a dense highly 
radio-absorbent regolith, and we vary the fraction of the part 
of the lunar surface relevant for lunar Cherenkov events due 
to neutrinos from Centaurus A and the Galactic Centre. For 
our "best-case regolith" estimate we use a 0% mare fraction, 
since there are very few mare regions in the lunar South 
where events from Centaurus A would be located. For our 
"worst-case regolith" estimate, we use a 30% mare fraction, 
which is approximately twice the lunar average and would 
correspond to the near side of the Moon being uniformly 
covered with mare. The mare regions in our "worst case re- 
golith" are modelled as having density p = 3.0 g/cro? and 
field attenuation length / = 2.94 m at 1 GHz, while "best 
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Figure 2. The exposure (km^-days) of the 2008 LUNASKA UHE u detection experiment usi ng the ATCA to 10^^ eV neutrinos. The 
small circles show the directions of UHE CR events above 5.6 X 10^^ eV detected by Auger ([Abraham et al.ll2007ah . The positions of 
the Galactic Centre and four important AGN are indicated by crosses labelled "Sgr A*", "Mrk 421", "Mrk 501", "M87" and "Cen A". 
The Galactic plane and Supergalactic plane are indicated by long-dashed and dotted curves, respectively. The region of sky accessible 
to lunar Cherenkov observations using the ATCA is between the two chain curves labelled "ATCA RANGE" . In the declination range 
— 10° < (5 <+15° (between the short dashed curves labelled "ANITA") the exposure of ANITA far exceeds that of other experiments, 
but for clarity has not been plotted. 



case regolith" is modelled with p — 1.8 and £ = 22 m. Fi- 
nally, we assume a (10 ± 10)% increase in aperture due to 
the interactions of seco ndary fi and r, as calculated for other 
similar experiments bv lJames fc Protherod l|2009ah . Our er- 
ror estimates do not account for the extreme uncertainty in 
the UHE neutrino interaction cross-sections - these are so 
poorly understood at these energies that determining them 
is an experimental goal as much as a calculation uncertainty. 
Given the strong proportionality between cross-section and 
effective area /sensitivity for this and similar experiments 
l|james fc Pro thcroc 2009a), our flux limit is really a limit 
on the flux — cross section product. 



4 SUMMARY AND CONCLUSION 

In Fig. Owe show the all-flavor neutrino flux limits for Cen- 
taurus A. With Centaurus A only 3.7 Mpc away, and with 
the pion photo-production energy-loss distance o n the CMB 
minim izing at ~12 Mpc above lO'^^ GeV (e.g. |Prothero3 
I2OO4I ') for rectilinear propagation one would observe UHE CR 
almost unatt enuated by pion photo-prod uction interactions 
on the CMB. lRieger fc Aharonicinl l|2009l ) suggest that shear 
acceleration along the kpc jet may accelerate protons be- 
yond 5 X 10^^ eV. The m agnetic field of Centaurus A's giant 
lobes l|Feain et al.ll2009l ) may also pr ovide an environment 
suitable for acceleration of U HE CR (|Benford fc Protherod 
I2OO8I: 'Hardcastle et al."200^). 

iCupco & Hanncstad (20ol) have predicted the flux of 
UHE neutrinos from the Centaurus A core ('CH08' in Fig. [3)) 
using a model of an optically thick p i on ph oto-production 
source described by iMannheim et al.l (|200ll ). They assume 
that accelerated cosmic ray protons are perfectly magnet- 
ically contained, and escape only through photo-hadronic 
interactions which convert them to neutrons. Under their 
model, the observed spectrum of cosmic rays requires 

a E~^''^ proton injection spectrum within the source, which 
would also produce a E~^'^ spectrum of neutrinos; however. 
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Figure 3. Neutrino flux limits for Centaurus A from our 2008 LU- 
NASKA ATCA experiments, and based on the GLUE, NuMoon 
and RICE, as well as published results from IceCube and Auger 
experiments. Limits are shown as integrated limits to a E^^ spec- 
trum (dashed lines), as well as model-independent limits (solid 
curves) where available - the NuMoon limit is model-independent. 
For the LUNASKA ATCA limit we show results from standard 
modeling (thick curve and thick dashed line) and using a toy 
model of small-scale lunar surface roughness (thi n curve and 
thin d ashed line) as described in Appendix B of I James et al.l 
||2010|) — the abrupt transition near 2 X 10^^ eV is a model arti- 
fact. Neutrino flux predictions of t wo AGN mod els for UHE CR 

roduction as labeled: KO T09 (.Kachelriess et al... 2009a') : CH08 

Cuoco fc Hannestadll200d ). 



these assumptions may break down towards the upper end 
of the energy range considered here. They normalize the CR 
flux by assuming that 2 of the Auger events above 5.6 x 10^^ 
eV are from Centaurus A, and determine the relative nor- 
malization between the CR and neutrino fluxes through 
Monte Carlo simulations of p-^ interactions in the accelera- 
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Table 1. Experimental exposures (km^-days) of GLUE, RICE, NuMoon (Westerbork) and the LUNASKA ATCA observations to UHE 
neutrinos at discrete energies from the Galactic Center and Ccntaurus A. 
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Figure 4. Neutrino flux limits for the Galactic Centre from our 
2008 LUNASKA ATCA experiments, and based on the GLUE, 
NuMoon and RICE experiments, and published IceCube limit. 
Neutrino limits are shown as integrated limits to a E^^ spectrum 
(dashed lines), as well as model-independent limits (solid curves) 
where available - the NuMoon limit is model-independent. Thin 
curves for the LUNASKA ATCA result are for modelling using 
a toy model of small-scale lunar roughness. SHDM models (see 
text) with Mx = lO" GeV and Mx = lO^'^ GeV (as labeled). 
Limit at 10^ GeV: Auger limit on the cosmic ray flux f rom a point 
source at the Galactic Centre ll Abraham et al.ll2007bl) . 



tion region using the SOPHIA event generator (|Miicke et aLl 
l2000f ). 

iKachelriess et all (|2009al ) also assume a CR flux com- 
posed of protons, normalized with the assumption that 2 of 
the UHE Auger events are from Centaurus A. They con- 
sider several possible proton injection spectra, with acceler- 
ation occurring either in regular electromagnetic fields close 
to the core of the AGN or through shock acceleration in 
the jets, and predict the resulting neutrino and gamma-ray 
spectra for each. The jet acceleration scenarios are excluded 
by TeV gamma-ray data (|Kachelriess et al]|2009bl ). and we 
plot ('KOT09' in Fig. [3]) their result for acceleration near 
the core with an proton injection spectrum. 

In Fig. 2] we show all-flavor neutrino flux limits for 
the Galactic Center. We consider the possibility that the 
Galajcy's dark matter halo is composed partly of super-heavy 
dark matter, which decays or annihilates into particles which 
cascade into neutrinos, photons, and nucleo ns. We use frag- 
mentation functions of lAloisio et all (|2004l ) for Mx = 10^" 



GeV and normalize the sum of the gamma-ray plus nucleon 
components to the Auger limit on t he cosmic ray flux from a 
point source at the Galactic Centre ([Abraham et al.ll2007bl ). 
We plot the corresponding neutrino flux shown as the curve 
labeled "10^'' GeV". We also show the expected neutrino 
flux assuming Mx = 10" GeV. Of course whether or not 
the neutrino (and cosmic ray) flux from SHDM annihilation 
appears point-like wil l depend on the r adial distribution of 
the dark matter halo (jEvans et al.P 2002). If this distribution 
is cusped the angular distribution will be narrow, approxi- 
mating to a point source, but in the case of other models the 
angular distribution can be broadened to up to ~60° half 
width. 

In comparison with experiments such as ANITA, our 
methods to improve sensitivity to certain speciflcally tar- 
geted regions were successful. We have reported the first 
neutrino fiux limits above 10^^ eV from Centaurus A and 
the Galactic Center. For Centaurus A, the limit from our 
2008 LUNASKA experiments using the ATCA is currently 
the strongest at primary neutrino energies of 5 x lO^^ eV 
and above. While not ruling out any of the current theo- 
ries of UHE CR production we have proved the viability of 
the lunar Cherenkov technique for targeted observation of 
specific UHE neutrino source candidates. 
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